Spin dynamics of two-dimensional electrons with Rashba spin-orbit coupling and 

electron-electron interactions 
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We study the spin dynamics of two dimensional electron gases (2DEGs) with Rashba spin-orbit 
coupling by taking account of electron-electron interactions. The diffusion equations for charge 
and spin densities are derived by making use of the path-integral approach and the quasiclassical 
Green's function. Analyzing the effect of the interactions, we show that the spin-relaxation time 
can be enhanced by the electron-electron interaction in the ballistic regime. 

PACS numbers: 72.25.Rb, 71.70.Ej, 71.10.Ca 



I. INTRODUCTION 

Spin based electronics [H or spintronics has been 
an active research area in the past decade. The effort 
for effectively manipulating electron spin by means of an 
applied electric field [1, [J, is an important issue there. 
The system with spin-orbit (SO) couplings makes those 
efforts possible and thus brings great interests from both 
academic and practical aspects recently. Thus, it is essen- 
tial to study the spin relaxation for further development 
of spintronics. 

There are four main mechanisms of spin relaxation in 
semiconductor systems 0, i, 0, i, i, EqI ■ In Elliott- Yafet 
mechanism, the spin-orbit coupling induces a mixing of 
wave functions for valence-band states and conduction- 
band states. The mixing that results in the spin relax- 
ation of electrons is due to the scattering by impurities or 
phonons. The Elliott- Yafet mechanism operates in semi- 
conductors with and without a center of inversion sym- 
metry, while it is most prominent in the centrosymmetric 
ones (such as silicon) . The Bir- Aronov-Pikus mechanism 
is applicable for p-doped semiconductors in which the 
electron spin flipping is induced by exchange interaction 
with holes. The hyperfine interaction provides another 
important mechanism for ensemble spin dephasing 
and single spin decoherence of localized electrons. The 
D'yakonov-Perel mechanism depicts that electrons can 
feel an effective random magnetic field arising from the 
spin-orbit coupling in systems with inversion asymme- 
try such that spin relaxation occurs. This mechanics 
can interpret the spin dephasing in crystals without in- 
version center and is particularly applicable for n-type 
samples. For two dimensional n-type semiconductor sys- 
tems without inversion symmetry, the D'yakonov-Perel 
mechanism is believed to be most important in wide 
ranges of carrier temperature and concentration. Under 
certain conditions, the Elliott- Yafet mechanism may af- 
fect spin dynamics of two-dimensional electrons in these 
systems. The Bir-Aronov-Pikus mechanism is impor- 
tant for p-type semiconductor systems and the hyperfine- 
interaction mechanism dominates for localized electrons. 

Most studies of spin relaxation in semiconductors have 



focused on impurity (somewhat less phonon) mediated 
spin flips while neglecting the effect of electron-electron 
interactions for a long time. It has been noticed re- 
cently that electron-electron interactions play certain role 
in spin relaxation and dephasing in semiconductor sys- 
tems. The electron-electron interaction is known to play 
a crucial role in determining the transport and thermo- 
dynamic properties near the metal-insulator transition in 
two-dimensional electron systems [12[ , which is suspected 
to affect the spin relaxation for the spin susceptibility 
behaving critically when the metal-insulator transition 
occur 

[H, Hi, El- There are several experimental and 
theoretical studies on the effect of electron-electron in- 
teractions on spin relaxation. The electron-electron scat- 
tering results in additional momentum relaxation which 
induces spin dephasing of electrons through the motional 
narrowing of the D'yakonov-Perel type [lal as measured 
in n-GaAs/AlGaAs quantum wells The electron- 

electron scattering effect on the spin dephasing has been 
considered[l^ in a magnetic field, and a momentum de- 
pendent effective random magnetic field induced by the 
electron-electron exchange interaction can lead to spin 
dephasing of electrons [lo, Hi], It is also observed 
that the spin relaxation caused by the D'yakonov-Perel 
mechanism gives considerably different rates depending 
on the technique employed [23|. 

However, as we are aware, the explicit form of the 
diffusion equation for two-dimensional electron gases 
(2DEGs) with spin-orbit couplings has not been derived 
by taking account of electron-electron interactions. It is 
thus obligatory to develop the explicit form of the dif- 
fusion equation to study the spin dynamics for 2DEGs 
with spin-orbit couplings as well as electron-electron in- 
teractions. In this paper, we focus attention on the 
D'yakonov-Perel spin-relaxation mechanism. We inves- 
tigate the spin dynamics of electrons in two-dimensional 
n-type semiconductor systems with electron-electron in- 
teractions and Rashba spin-orbit coupling. 

The paper is organized as follows. In Sec. [Hi we 
take account of the electron-electron interaction for the 
2DEGs with the Rashba spin-orbit coupling. Applying 
the path integral formulation, we decoupled the interac- 
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FIG. 1: The Keldysh contour. 



tion in terms of an auxiliary Bose field. In Sec. IIIIl we 
employ the quasiclassical Green's function to investigate 
the spin dynamics of electrons. In Sec. IIV( the diffusion 
equations for spin and charge densities as well as the ex- 
plicit expression of spin-relaxation time are derived. A 
summary is given in Sec. |V] and some complicated for- 
mulae are given in the Appendix. 



II. AUXILIARY FIELDS DESCRIBING THE 
ELECTRON-ELECTRON INTERACTION 

Taking the electron-electron interaction into account, 
we study the spin dynamics of electrons in two- 
dimensional systems with structure inversion asymmetry. 
As the Fourier transform of the Coulomb repulsion be- 
tween electrons reads V{q) = 27re^/|q|, the Hamiltonian 
of such a system is given by 



H 



A, A' 



1 

+ 1 



(1) 



where -0'^^ (r) and -0'^ (r) represent the field operators with 
A =t,i labelling the spin state of the electron, /5(q) 
represents the Fourier transform of the density operator 

p(r) = X^A and a = {<7x, (Jy, (Jz) the Pauh 

matrices in spin space, U{r) a random disorder poten- 
tial and n the chemical potential. The other notions in 
Eq. ^ are A — with L referring to the size of the 
sample and b = ap x with a referring to the Rashba 
spin-orbit coupling strength. In holonomy representation 
(or the called coherent state representation), the Green's 
function can be expressed as a functional integral over 
the Grassmann fields and that reflect the fermionic 
nature of electrons. 



Here we adopted the unit h = 1 and the simplified no- 
tation Dtp = Dtp^ DipK The action S'fV', V'] in the above 
equation is given by 



dt[ I dV^^^(r,t)VFAA'V'^'(r,i) 

AA' 



where P^AA' = {-id / dt-Vl/2m+U {r)- p)5xx' +h-axx- ■ 

We divide the fermionic field t) into two compo- 

nents V'i(r, i) and -02(r, i) which reside, respectively, on 
the upper and lower branches of the Keldysh time con- 
tour shown in Fig. ([T|). Hence, the second Hue of Eq. ([3]), 
which refers to the interaction part, can be written as 

<S'mt[V'l,'0l] - Sint[^2,'4>2] with 

/2 



where i = 1,2. With the help of two auxiliary bosonic 
fields ^i(r,t), we can decouple those two terms rele- 
vant to electron-electron interactions via the Hubbard- 
Stratonovich transformation [24] , namely, 



/2 

J D4>,{q,t) exp[^ J di^^^.(q, t)M~^, t)] 



X exp 



q#0 



+p,{q, t)4>^{-ci, t)} 



Then we can write the Green's function as follows, 

Gxx'{r,t;r',t') 

J D^D^D^'f^ir, t)^^' (/, t') e-*^[*'*'*l 



(4) 



(2) where the action in real space is given by 
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XX' 

+ j dtj dW/{<i>^(r,i)^Fo-i(r-r')fT3$(r',i)}, 



(5) 



r 



in which the Pauli matrix (73 = diag(l, —1) is defined on 
the Keldysh space, and Vj^"^ is defined via the following 
relation 

Jd^iVoir ~ ri) - r') = S{r - r'). 

The other notions appeared in Eq. (O are fermionic dou- 
blet V?, bosonic doublet $, and vertex matrices 7^, they 
are defined as 



^1' 
^2' 



1 \ , / 

^1 = I ' ^2= 0-1 



For calculation convenience, one can introduce a parti- 
tion function for the coupling between the fermionic and 
bosonic doublets, 

= ( Tc e-''^«[*'*l 

A 

where Tc stands for time ordering along the contour C 
and ( • • • )>f means functional integration over ^E* field 
with the action 

5[*]= / dtd2r{^TyAA'*^r,i)a3*^'(r,i)}. 

AA' 

Then the Green's function in Eq. ^ can be formally 
expressed as a functional integration over the bosonic 
fields. 



Gxx'{r,t;r',t') = Nj D<i> Gxy{r,t;r' ,t' \ $) 

xexp{ (6) 

where the normalization coefficient is denoted by N and 
the action Se[^] is defined by 

Se[<i>] =i\nZ[<i>] + j dt dPrdPr' 

x{$^(r,t)^V'(r-r>3'i>(r',t)}, (7) 

and the kernel (^(r, t\ r', t' \ $) is given by 
Gxx'{v,t-v',t' I $) 



We can average the Green's function Gw (r, t; r', t') over 
disorder as follows (25j 

{Gxx'{v,t■,v\t'))d^s = N J D<i>{Gxx'ir,t;r',t' \^))d^s 

xexp{ (9) 

where (• • • )dis refers to the average over disorder. The 
random disorder potential U{r) is assumed to be charac- 
terized by a correlation function 

1 



{U{r)U{r'))d^s 



2'KVT 



where v — m/Trh^ stands for the density of states. The 
average of the Green's function over disorder introduces 
the elastic scattering time r which is relevant to the ran- 
dom disorder. We neglect correlations between the meso- 
scopic fluctuations of {Se[^])dis and the fermionic opera- 
tors in Eq. ^ so that the average of the Green's function 
G\y (r, t; r' ,t' \ <i>) can be separated from the bosonic ac- 
tion 5'e[^']. This approximation is valid since the meso- 
scopic fluctuation is smaller than average quantities. 

After averaging over disorder, we rotate the Keldysh 
bases G — > La^GL'^ through a unitary matrix L 

so that the Green's function takes the following shape, 
G(r,i;r',t'|$) 

_ f G"{r,t-Y',t' I G^{v,t;Y',t' \ $ 



L 



\ G^{r,t;r',t' \ $) G^{r,t;r\t' \ $) 



(10) 



Note that the Green's function G(r, t; r', | $) is a 2 x 2 
matrix defined in the Keldysh space, of which the matrix 
entities are again 2x2 matrices defined in spin space. 
The bosonic fields after rotation take the following two 

components (pi — ^('^'^ '^2) ^^'^ = 2^'^^ ~ "^2) that 
reside on the upper and lower branches of the contour 
C, respectively. Then the corresponding vertex matrices 
turn to 7i(2) = L{^i ± 72)o-3Lt, namely. 



71 



and the interaction term ^''^(—60^7^)^'^ becomes 















, 72= (5 





As 71 and 72 constitute a repre- 



sentation of Z2 group, the interaction can be regarded as 
the coupling between Fermi field and Z2 Bose field. 
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These Bose fields define the following propagators: 



equations in the saddle point approximation, namely, 



D^{ri,r2;h,t2) = -2^(0l(rl, ii)0i(r2, ta)), 
Z?^(ri,r2;ti,t2) = -2i(0i(ri, ti)02(r2, ^2)), 
I?-^(ri,r2;ti,i2) = -2i{M^i^h)M^2,t2)), 
(02(ri,ii)02(r2,i2)) =0. (11) 

One can show that those propagators obey the Dyson 



(12) 



with notations x = (r,t), -D(f(q) = DQ{q) 
and 



-2^eV|q|, 



J 



D = 



D"^ 



^" - ' 



n 







n^(xi,a;2) = n^(a;2,a;i) 
n^(a:i,a;2) 



JTrsG^(ri;ti,ti | $) 
-2z ,50i(r2,i2) ' 
(5Tr,[G^(ri;ti,ti | $) + G^(ri; ti, ii | $)] 



-2z (502(r2,t2) 



(13) 



where Trg stands for the trace in the spin space. 



III. THE KINETIC EQUATION 

In previous section, the electron-electron interaction 
has been decoupled with the help of auxiliary bosonic 
fields 4'i{2)- This means that the influence of the interac- 
tion can be described by a Z2 Bose field, 



4>i{r,t) 02(r,t) 
02(r,i) 0i(r,t) 



Now we are able to apply the quasiclassical Green's func- 
tion ^,t23i[2l| approach to study the spin dynamics. We 
derive the Eilenberger equation from the right-hand and 
left-hand Dyson equations obeyed by the Green's func- 
tion G(r,t;r',i' | $) in Eq. HI 



dtg + VF • V5 + z [b • (7, g] = 



.9 { 9 



9 )n 9 



2 T 



(14) 



where vp denotes the Fermi velocity, r is the elastic scat- 
tering time arising from the adoption of the standard 
self-consistent Born approximation; (•••)n means tak- 
ing average over the direction of the electron momentum 
^ = P/|p| = {cos9, sin 6*), and the covariant derivative 
is defined by 

dtg = dt,g + dtj + iipir,ti)g - igip{r,t2). (15) 

The quasiclassical Green's function in Keldysh and spin 
spaces. 



9' 9^ 



(16) 



can be derived by integrating the Fourier transform of 
the Green's function in Eq. (jlOp over energy variables, 
i.e.. 



){ti,t2;n,v) = - d(,G{ti,t2;p,r), 



G{h,t2;p,r) = y dV'e^P'-'G(ri,ii;r2,t2 | $), (17) 



where ^ = p^/2m — /i, r' = ri — r2, r = (ri + r2)/2. The 
electron polarization operators can be obtained in terms 
of Eq. ^ and Eq. i.e.. 



n^(a;i,a;2) = n^(a;2,a;i) 

fde. TT STY,g^{ti,h;n,r^) ^ 

-J + 2 5Mr.,t2) ]' 

n^(a;i,X2) 

deSTT,[g'^{h,h;n,ri)+g^{h,ti; n, n ) 
2n 2<502(r2,t2) 



(18) 



Since physical observables are determined by the 
Keldysh component of the quasiclassical Green's func- 
tion, namely ( g^ {ti,t2;n,r) )$ (here the subscript $ 
refers that the functional average [2^ is taken over the 
field $), we need to solve this component from the 
Eilenberger equation. Decomposing the Green's func- 
tion in charge and spin components ( g^ (ti, t2; n, r) )$ = 
9o + S^' ' '^1 one can obtain the charge and spin densities. 
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respectively, 

S(r,i) = -ii. J de(g^(t,e;n,r))„. 



where the diffusion propagator Tp is defined by 



(19) 



Now we turn to the kinetic equations for the two inde- 
pendent components and . For ( t;'^ )$ = in all 
orders of the perturbation theory, we have 

where the fluctuation parts 6g imply the effects con- 
tributed by the auxiliary bosonic fields. One can obtain 
from Eq. (|14p that Sg^ obeys the following equation, 

{dt + • V)<5.g^ + i[h-<j, 5g^] - - [5g^ - {5g^)^] 



-2i02(r,ti)^(ti -t2)Is, 



(20) 



where Ig denotes the unit matrix in spin space. When 
deriving the above equation, we have used the conditions 
gR = 6{h - t2)Is - g''Sg''/2 and g^ = ^6{h - t2)Is + 
5g^g^ /2. Equation gives rise to 

/f do' 
dridh / ^ 

X02(ri,i3)rp(t3 - ti,n',n; ri,r), 
r,(i,„',„;ri,r2)=./^e-(--r.)-..* 

xrp(n',n;c^,q), (21) 



{-iuj + ivFn- q)Tp{n,n';uj,q) + - [Tp{n, n';a;,q) 
-(rp(n,n';cc>,q))„] =2^<5(n-n'). (22) 



After obtaining the explicit form of Sg^ , we can further 
solve the 6g^ from the following relation 



(dt + VF • V)<5.g^ +i[h-<j, + i [Sg"" - ((5g^)„] 
= 2iMr, h)5{ti - t2)Is - « [01 (r, h) - 0i(r, ia)] (.9^)$ 



We take only the zeroth and first angular harmonics 
into account in the Keldysh component assumed spatial 
smoothness, 

( g^(ti,i2;n,r) )$ w ( g^(ii,t2;n',r) 

+2n- ( n'5^(ti,t2;n',r) )$^„,. (24) 

Decomposing the fluctuating term in charge and spin 
components 5g^ — SgQ + Sg^ ■ a, one can easily ob- 
tain the explicit expression of the 6g^ which is given in 
Eq. (jAl|l . The fluctuation part 5g^ related to the spin 
components fulfils the following equation. 



(9t -f VF • V),5g^' - 2b X + ^[^g^ - ((Jg^On] = -«[0i(r,ti) - 0i(r,t2)] ((g^(ii, t2; n, r))*,„ 
+2n. (n'g^(ti,<2;n',r))*,„,) + ±{{g^)^^^{{6g^)^ - Sg^){g'')^,n + (g^)*,„((<5g^)„ - 6g^) {g^)^,^} . (25) 



If denoting 




we can write Eq. (|25p in the following matrix equation, 
{dt + VF • V)Q - 2C0 + - (Q - (Q)n) 

T 

= -i[0i(r, ti) - 0i(r, f2)] ((-^fe)*,n + 2n • (n'Lfe)$,n') 

+ (^fc)*,n((5ff^)„-<5ff^)(5o''>*,n}, (26) 



where the matrix C is given by 

(0 —app cos6 \ 

-apF sine . 

app cos 9 appsinO / 

Then equation (|26p can be solved by utilizing the follow- 
ing expression 

{-iu! + ivpn • q - 2C)rs(n, n'; q) + -[rs(n, n';uj,q) 

-(r,(n,n';c^,q))„] = 27r(5(n - n'). (27) 

We give the explicit expression of Sg^ in Eq. (jA2p in the 
appendix |^ 

Since the concrete forms of 6g^ and 6g^ have been 
obtained, we can write down the kinetic equation satis- 
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fied by the Keldysh function through averaging the K 
component of Eq. over the auxihary bosonic fields, 



and the elastic collision integral is given by 



= Ce/{(.g^)*}+a„{(5^')$}, (28) 
where the inelastic collision integral reads 
Cm{(.g^)$}(ii,t2;n,r) 

= -z( [(/.i(r,ti)-0i(r,t2)]<5g^ (29) 



where 



Cei{{9'').s>}{tut2;n,r) = i[( (^^(ti, ia; n, r))* )„ - (.g^'(ti, tz; n, r))$] 

/dO 
dt3^[{g'^{ti,t3; ni, r))<E,A^(t3, ^2; ni, n; r) - (.g^(ti, ts; n, r))<E,A'^(i3, ^2; n, m; r)] 

d<3^ [A^(ii, ts; n, m; r)((7^(i3, ^2; ni, r))$ - A^(ti, ia; m, n; r)((7^(t3, i2; n, r))$] , (30) 

A^(ti,t2;n,ni;r) = ^ J dtsi [(5g^(ti, t3; m, r) - ^^^(ti, ^3; n, r)] 5g^(i3, tz; n, r) )<i., 
A^(ii,t2;n,ni;r) = ^ J dhiSg'^ {ti,h-n,r)[5g'' {t3,t2;ni,r) ~ 6g^ {t3,t2;n,r)])^. (31) 



One can see from Eq. ([29]) that the right-hand side is 
always zero. Thus we obtain 



r 

Substituting the explicit forms of 5g^ and Sg^ given in 
Eq. dni), (ED) and (|X2|) into Eq. one get the ki- 

netic equation which can be used to study the influence 
of electron-electron interaction on the spin dynamics of 
2DEGs with Rashba spin-orbit coupling. After some te- 
dious calculation, we obtain the explicit expressions of ^j^^^ -^^^-^^ ^^^^ ^^^^ ^^^^1 ^^^^^^ ^^^^^ 
the inelastic and elastic collision integrals, respectively, 
which are given in Eq. (|A3p and (|A6p in the appendix [XI 



deC™{(.g^)*}(t,e;n,r) =0. 



(33) 



IV. SPIN DYNAMICS 

After taking average over the direction of the momen- 
tum n, one can see from Eq. (|29II30|) that the elastic col- 
lision integral vanishes. 



is conserved, but also the number of electrons moving 
along a concrete direction n is conserved. 

Decomposing the Green's function in charge and spin 
components in the approximation of Eq. ([M]) . separat- 
ing the zeroth and first angular harmonics and utilizing 
Eq. ([32]) and Eq. we obtain from Eq. ^ that 



WFn- V(5^(t, e;n,r))4 



(32) 



{g''{t,e;n,r) = Cei{(5^)$} (i, e; n, r), (34) 

dt{g^{t, e; n, r))$,„ +i[h-a,2n- {n'g^{t, e; n', r))$,„,] 

, . , 1,. . . , 1 -K«Fn.V(2n.(n'.g^(t,e;n',r) =0. (35) 

but the average 01 the melastic collision integral over the 

direction does not vanish. This means that the elastic There is no contribution of the inelastic collision integral 

collision integral preserves the number of electrons on a to the spin dynamics due to the condition Eq. Solv- 



/ ^^ei{(.9^)*}(i,e;n,r)=0. 



given energy sheU defined by Eq. dM]), while the inelastic ing Eq. ((Ml) and substituting {n'g^(t, e; n', r))<E,,„' into 
collision integral does not preserve it. When ti — t2, 
equation (|A4[) gives rise to 



j deCin{{g'^)i>]{ti,c,n,v) = C„{(g^)$}(ti, ii; n, r). 



Eq. (j35p . we can obtain the spin and charge diffusion 
equations. The diffusion equation for the charge density 
reads. 



(36) 
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where dx = {dx,dy) and Co = v'ft/2. We introduce 
the distribution function / which reduces to the Fermi 
distribution in equiUbrium, 



/ = /o + a-fk = i(l-i/^'). 



(37) 



In the time r, the charge density becomes isotropic but 
the spin relaxation process does not start yet, hence [sO] 



given by 



M 



X Im(i?r)5(f (e) + - f-ie))Rr9^ie - to)] }, 

X Im(i?f )5o^(e) + (Uie) - /_(e))i?f ^^''(^ " ^)] }, 



1 1 



go^(e) = 2(l-2/o(e)), 
g^^(t,e;r) --4fk(t,e;r), 



(38) 



where 



/o(6) = (/+(e) + /_(6))/2, 
fk(i,e;r) = (/+(e)-/_(e))s(t,r), 

/,(e) = [exp(i^) + r\ 



(39) 



where s = {sxt Sy, Sz) denotes the unit vector along the 
spin, /±(e) represent the distribution functions projected 
along the direction parallel or antiparallel to the unit vec- 
tor s (all the energies are counted from the Fermi energy) 
and A/i = (/i-|_ — /i_) refers to the difference between the 
chemical potentials /i± of the electron spin subsystems. 
The diffusion equations for the spin components are given 
by 



dtSx — Cnd^Sx — 2CEdxSz H — jSx 

n 



~~~Sx + Fx{Sx, Sy, Sz), 



dtSy — Cod^Sy — 2CEdySz H — jSy 

— ~~^^y + Fy{Sx, Sy, Sz), 

T 

yy 

2 2 
dtSz — Cody^Sz + 2CEdxSx + 2CEdySy H — ^5*2 

1 



~~r^z Fz{Sxt Sy, Sz), 



(40) 



2{apF?T 
M 



X (lm(i?r)+Im(i?|^))5o^'(6) 

+ (/+(£) - f-{e)) [RT + RT)9n^ - ^)] }, 



(41) 



where M — J (ie(/+(e) — /-(e))- In order to obtain the 

concrete expressions of the characteristic times t^^, we 
firstly take the energy integration in Eq. (|¥l]) . Since the 
spin splitting is small, i.e.. 



fi+ - i-i- |<| fi+ I, I I, 



the energy integration can be taken as follows. 



M 



de{f+{e -uj)-Ue~ oj)) (/+(e) + /_(e)) 



2 1 de'-Mp^m 

oe 



de 



dfoje) 
de 



9 , , ^ N 

1 — — — (wcoth— — — ). 
du^ 2kBT' 



(42) 



After the energy integration, the characteristic times t| 
have the forms 



Siappfr I — [_(^coth— -)] 



TSx Tyy ' ' Jo ^TT ^ duj 2kBT 



[im(i?r)-i?r], 



(43) 



where Ce = avpPFT, = l/[2(apF)^T] and 

Fi{Sx, Sy, Sz) is a quadratic form of {Sx, Sy, Sz) lack- 
ing of the Sj {£ = X, y, z) term. The characteristic times 
rff^ describe the effect of the electron-electron interaction 
on the spin relaxation, and their explicit expressions are 



the detail of the calculations of the kernels FC^ and Imi?2"' 
are given in appendix B. Now we discuss the influence of 
the electron-electron interaction on the spin-relaxation 
time in the ballistic regime Tr ^ 1. 

We can obtain the characteristic time in the ballis- 
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tic regime utilizing the kernels and Imi?2"' in Eq- (|B7p 

-(Tr»l) = 8(.p.)Vy^ -[-(.coth^)] 
I- —1 .Stt _-, 2UJT 

= ^(^^(2fc.T-i^.coth-^), (44) 

where tan"""^ (ojt) is replaced by 7r/2 for ujt ^ 1 in the bal- 
listic regime and Ep is in the place of the upper limit of 
the integral. In the low temperature regime fcsT <C Ep, 
the second term approaches a constant independent of 
the temperature, so the first term manifests the temper- 
ature effect in the contribution of the electron-electron 
interaction to the spin-relaxation time. 

When the total spin density S is spatially homo- 
geneous and parallel to the £th-axis of the coordi- 
nate frame, the contribution of Fi{Sx, Sy, Sz) vanishes, 
namely Fi{Sx, Sy, Sz) = 0. The diffusion equations for 
spin components Si can be simplified, for example, 

dtSx^~^Sx + ^Sx^-^Sx, (45) 

' s ' XX ' XX 

t' 

where = t^/(1 —). Therefore, the spin-relaxation 

XX 

times can be determined by and r|^, consequently. 



' XX ~-i ^ yy ' 

1- ^ 1- ^ 

XX yy 

(.rLY' = {rLr' + {rly)-\ (46) 

We can see that the total spin decays exponentially when 
< T'^jrlg < 1. In terms of the explicit forms of 
the characteristic times t^^ in the ballistic regime, the 
spin-relaxation times involving the effect of the electron- 
electron interaction take the following forms 

t' 

where Tp = Ep/ks is the Fermi temperature. It is 
worthwhile to indicate that there exists an obvious en- 
hancement of the spin-relaxation time with increment of 
the temperature in the ballistic regime. The increasing 



amplitude of the spin-relaxation time depends on the ra- 
tio of the temperature to the Fermi temperature. In con- 
clusion, an obvious enhancement of the spin-relaxation 
time can be induced by the electron-electron interaction 
in the ballistic regime for systems under consideration. 

V. SUMMARY 

In the above, we presented a theoretical study of the 
influence of electron-electron interactions on the spin dy- 
namics for 2DEGs with Rashba spin-orbit coupling. We 
employed the path-integral approach and the quasiclas- 
sical Green's function to deal with the electron-electron 
interaction. With the help of the auxiliary Bose field, 
the electron-electron interaction was decoupled via the 
Hubbard- Stratonovich transformation. Then one is able 
to derive the Elienberger equation by using the Green's 
function after the transformation. Through tedious cal- 
culation, wc further derived the spin and charge diffu- 
sion equations, from which the spin-relaxation time can 
be given explicitly. We analyzed the influence of the 
electron-electron interaction on the spin-relaxation time 
in the ballistic regime and found an obvious enhance- 
ment of the spin-relaxation time with the increment of 
the temperature T. The increasing amplitude of the 
spin-relaxation time depends on the ratio of the temper- 
ature to the Fermi temperature. The electron-electron 
interaction changes the wave vector k and hence results 
in the variation of the spin precession vector. This ex- 
hibits that the electron-electron interaction plays an im- 
portant role in the spin relaxation of electrons when 
the D'yakonov-Perel spin relaxation mechanism domi- 
nates. It is expected to be helpful for understanding 
the spin dynamics of 2DEGs with spin-orbit couplings 
and electron-electron interactions. Our formulation can 
also be extend to the case of bulk inversion asymme- 
try, namely the additional Dresselhaus term i3l| with 
b = P{px,-Py) + j{PxpI - PyPD- 
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APPENDIX A: EXPLICIT FORMS 

The explicit form of Sg^ is 
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dt0[(t)i{Yi,ti-'te)~ (j)i{Yi,t2-te)\ I —Tp{te,n,YV]Y,Yi) 

/dO' dO" f ( 
2^ 27 / ^^'^i'^^9rp(ie,n,n';r,ri)|2i(/.2(ri,ii -te)J(ti -iz) 

+ ^[{ ^p{ti - i3,n",ni;r2,ri) )„i - T p{ti ~ is, n", n; r2, ri)] 02(r2, ^4) 
+ ( gxih - te,i3;ni,r) • ( g;^(i3,i2 - ie;ni,r) )<i>,„J j. (Al) 

The explicit expression of (5g^^ is 

/f do' 
dte[(j)i{vi,ti-te)-(l)i{vi,t2~t0)\ j — rs(ie, n, n'; r, n) 

X |( ife(ii - iff, t2 - ie; ni, r) )$,„i + 2n' • (niLfe(ii - te, i2 - ie; ni, r))$,„j | 

2i /■ dd' dd" f ,2 , ^ , , ^, 

+ ~ / '2^^ d ridteTs{tg,n,n;r,ri){ (ti - t^, ^3; rii, r) )$,„^ 

X { [( ^^(^4 - i3,n",ni;r2,ri) )„i - Tp{t4 - t3, n", ni; r2, ri)] 

X (I32ir4,t2){ Lk{t3,t2 ~te;ni,r) )^^-n^y (A2) 

The inelastic collision integral reads 
a„{(5'^)$}(t,e;n,r) 

X [( rp(w; ra, ri) ) ( T p{~uj; r, ri) ) - ( T p{-uj] r, ri)rp(tj; ra, ri) )] 

X [( g^it-, e - n, r))<i>^„ {g^ {t, e; n, r) )*,„ + ( g^(t, e - n, r) )<i.^„ • ( g^(t, e; n, r) 



y d^i y ^L»^(w;r,ri)cr™[( r,(-w;r,ri) ) + ( r,,(w;ri,r) )] [( e; n, r) )$,„-( Lfe(i, e - w; n, r) 
+ 1 1 d\,d\2j '^[D^{u:-v,V2)^D^{uj-V2,v)]ara 

X [(r,(-tj; ra, ri)) ( T p(-uj- r, n)) - (r,(-tj; r, ri)rp(tj; ra, ri))] (.g^(t, e - w; n, r))<i,,„(Lfc(i, e; n, r))<i>,„, (A3) 



where 

Cm{(g^)$}(ii,t2;n,r) 
= / ^C,„{(.9^)*}(^,.;n,r)e-(*-'^),(A4) 

and (Tp means angular averaging defined in 
Eq. (jAlOp . the matrix ct„i — {<Jx,o-y,o'z) and the tem- 
poral transformation of the Green's function has been 



used due to a much faster dependence on the difference 
ti — t2 than on the ti + ^2 



g^it^,t2; n, r) = I ^g^('^,e; n, r)e"(*-*^)(A5) 

the propagators of auxiliary fields have the same trans- 
formation. The elastic collision integral can be written 
as 
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Cei{(ff^)<i.}(i,e;n,r) = -^n,( n.g'^ {t,e;n,r) )<[,,„ - ^ J ^n,R\^{uj){ g,f (t, e - w; n, r) 



n^g (t,e;n,r) 

— ] 

27r 



do; 
2^ 



n,i?2^(cj)( nj5r^(t,e- w;n,r) g^(t,e;n,r) 



n^{Rl\uj)y{ g^(i,e;n,r) njg^(t, e - n, r) 

cr,„i?3^(t^)( ifc(t, e-cj;n,r) nj5^(t, e; n, r) 

nj5^(t, e;n, r) cr„i(i?3^(cj))*( Lfc(t, e - w; n, r) ■ 

I 



H — / — n 

T. 2tt 



(A6) 



Where n^j^) refers to the coraponent of the unit vector with (j)g being the angle between the wave vector q and 
n, with t,j — x,y and the kernels R\^(oj) — R^^it^) in the x-axis. In terms of the explicit form of the diffusion 

propagator, one can obtain 



Eq. (|A6P are defined by 



^Im J ^D^{u;;<D{{ rp(n,c.;q)n, ) 

x( n,r, ) - ^s,.,{{ rp){rp ) - ( r,r, ))}, (A?) 



rpn,r,n, )-(rpn, )(rpn, )}, (A8) 



1 



i?^(cc>;q){(r,n, )(n,r, ) 



-<5,,,((r, )(r, )-(r,r, ))}, 

where we have introduced the notation 
dOdO' 



(A9) 



(27r)^ 



/(n)r,(,)(n,n';o.,q)Mn')-(A10) 



APPENDIX B: CALCULATION OF THE 
KERNELS W 



Tp in Eq. ((22)) . we can obtain 

rp(n,n';a;,q) = 27r(5(n - n')rpo(n; w, q) 
+ rpo(n;w,q)rpo(n';cj,q) 



(Bl) 



where 



rpo(n;cj,q) = 



1 



-iu + ivpTo. • q + 7 
1 

-iuj + ivpq cos{4> — 



{-iuj + -)2 + vjpq'^, 

T 



(B2) 



(r,> 

(Epriy) 
(rpFp) 
(rpn^rEp) 



Tt- 1' 



^'""^"^ " ivpgiTr - 1) 



{nyTp) 



(T + iu) ) cos 4>„, 

T 

T 1 

— -(T + iuj- -) sin0g, 

ii'Fqy 1 T — 1) r 



T 



T 



i;^g2(Tr-l) 



(1 



(1- 



^ ^ )^cos', 



-^)2sin2 



According to the definition of the diffusion propagator p^x^ pnx 



Tt 



-ILO +7.2 

^ sm'' ( 



-^2 



(rpnyTpTij^ 



T-(^^+i) 
T2(Tr- 1) 

,-iuj + i „ 
■ cos 



■ COS 



Tt - 1 ^ T2 



T-(^^+i) . 
T2(Tt- 1) 



(B3) 



Utilizing above formulas, we find that the kernels RY{(^) 
and i?2"'(w) are diagonal, i?*"' = S^jRi, which can be writ- 
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ten as 
R 



Imi?2H = Im r ^D''iuj;q)\^ 
Jo 47r I^Wpg^ 



T(T- i)2 i' 
1 [T + zu;-il2 



g2 T(T-i) 
+ T(T_1)2 }- (B4) 

It is not difficult to calculate the concrete forms of the 
electron polarization operators from Eq. (jlSp . for exam- 
ple, 



^ _1_ -I 



(B5) 



Substituting the polarization operator into Eq. p2l) . we 



obtain the propagator of the Bose fields, i.e., 



1 r 



1 



(B6) 



where the approximation in the second line corresponds 
to the unitary limit associating with larger distances than 
the screening radius. 

We obtain the concrete expressions of the kernels Ri 
and i?2 in the ballistic regime Tr ^ 1, 

Ri{Tt > 1) (X 2"("T +tan"-^wT 



1+Cj2t2'' 



Imj?2(TT > 1) (X - + tan ^ - 



(B7) 
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